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We calculate holographically arbitrary n-point correlators of the boundary stress tensor in three- 
dimensional Einstein gravity with negative or vanishing cosmological constant. We provide explicit 
expressions up to 5-point (connected) correlators and show consistency with the Galilean conformal 
field theory Ward identities and recursion relations of correlators, which we derive. This provides a 
novel check of flat space holography in three dimensions. 
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INTRODUCTION 


Correlation functions are the main observational and 
theoretical entities in many branches of physics, in par¬ 
ticular in any system that effectively is described by some 
quantum field theory (QFT), with numerous applica¬ 
tions in condensed matter physics, cosmology, particle 
physics and many other research fields. Given their ubi¬ 
quitousness, it is useful to have several tools available to 
calculate correlators. 

A remarkable tool for strongly coupled QFTs is pro¬ 
vided by holography [1, 2]. In particular, the Anti-de Sit¬ 
ter/conformal field theory (AdS/CFT) correspondence 
[3, 4] allows to calculate correlators in a specific strongly 
coupled CFT by mapping this (typically) hard calcula¬ 
tion to a simple one in a specific gravity theory [5, 6]. 

As the simplest example, the 1-point function [i.e., the 
vacuum expectation value (vev)] of the stress tensor T ^ 
is calculated by taking the first variation of the corre¬ 
sponding gravity action S Ads with respect to the metric, 
evaluated on solutions of the equations of motion (EOM): 
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Note that the right hand side is non-zero and finite due 
to boundary terms [7-9]. The holographic calculation 
of arbitrary (connected) n-point correlators of the stress 
tensor is a straightforward generalization of (1) and re¬ 
quires to evaluate the n th variation of the gravity action. 
The n-point correlators of the stress tensor can also be 
calculated on the CFT side, so ultimately these calcula¬ 
tions are a check of the holographic Ward identities. 

For more recent holographic correspondences that go 
beyond AdS/CFT, such as flat space holography [see [10- 
34] for some recent advances in three dimensions (3D)], 
similar checks would contribute to increase the credibil¬ 
ity of the purported relations between certain QFTs and 
their conjectured gravity duals. 

The main goal of the present work is to provide the 
check outlined above for flat space Einstein-Hilbert (EH) 


gravity in 3D by calculating arbitrary n-point correlators 
of the stress tensor components on the QFT side and, 
independently, on the gravity side. So far only the 0- 
and 1-point functions were calculated holographically [18, 
21, 22, 25]. Our results are the first explicit check of 
the flat space holographic 2- and higher-point correlators 
and allow to present for the first time explicit results 
for the 3-, 4- and 5-point correlators of the stress tensor 
components in two-dimensional (2D) Galilean conformal 
field theories (GCFTs) as well as their Ward identities. 


ADS3/CFT2 

We start with the AdSs/CFT 2 derivation of n-point 
correlators to exhibit the required tools in a familiar con¬ 
text. The key result is that all n-point correlators of the 
stress tensor in a CFT 2 are determined by the underly¬ 
ing symmetry algebra, which consists of two copies of the 
Virasoro algebra (n, m are integers), 

[Gi ^m] = (ti rr/jCn+m -\- :~n(n l)J n _ m (2) 

and similarly for the barred generators C n , with the cen¬ 
tral charge c replaced by its barred counterpart c. (The 
same remark applies to all remaining CFT formulas, so 
to reduce clutter we display only unbarred quantities.) 

In their seminal paper, Belavin, Polyakov and 
Zamolodchikov proved a recursion relation between n- 
and (n — l)-point correlators of the stress tensor in any 
CFT 2 [35]. For the connected part of the correlators this 
recursion relation reads 

(T 1 T 2 ... T") = V ( 4- + — d A {T 2 ... T n ) (3) 

7^2 Wii / 

where we introduced the abbreviation zij = Zi— Zj. Here 
T l = T zz (zi) is the holomorphic component of the stress 
tensor and z, z are coordinates on the plane with metric 
9zz = 0 = g zz and g zz = 1. The disconnected contri¬ 
bution to the n-point correlators is the sum of all per¬ 
mutations of products of lower-point functions such that 
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all products have n points. Since holographic computa¬ 
tions yield connected correlators we focus on them. All 
correlators in this paper are solely the connected part. 

The n-point correlators are completely determined by 
the Ward identities and the 2-point correlators, which 
in turn are uniquely determined by conformal invariance 
( 2 ): 


(T 1 T 2 ) = 


9 r 4 
ZZ 12 


(4) 


Using this in the recursion relation (3) gives the higher- 
point correlators [36-39]. 
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This procedure can be iterated to arbitrary n-point 
correlators, with increasingly lengthy expressions. Thus, 
if we can holographically establish the recursion relation 
(3) and calculate the 2-point correlators (4) we have es¬ 
sentially succeeded in holographically computing all the 
n-point correlators of the stress tensor. To this end, 
we consider a deformation of a free CFT action Sq by 
a source term p for the stress tensor, 


Sfi = S 0 - J d 2 zp(z , z)T(z). 


( 6 ) 


We then localize the source at some point (z 2 , -22)- 

p(z,z) = e 6 ^ 2 \z - z 2 ,z - z 2 ) (7) 

Here e is a small expansion parameter introduced to keep 
track of the order of /j. The 1-point function with respect 
to the //-deformed vacuum yields the 2 -point correlator 
with respect to the free vacuum, 

(T 1 )^ = (T 1 ) + e (T 1 T 2 ) + 0(e 2 ) . ( 8 ) 

We now compute this expansion in the dual AdS 3 grav¬ 
ity theory. According to the AdS/CFT dictionary, the 


source p corresponds to the boundary value of the met¬ 
ric component g 2 ~. To compute the n-point function we 
need the n th variation of the bulk action with respect to 
the source. Such a computation is lengthy in the metric 
formulation, but fortunately simplifies significantly in the 
Chern -Simons (CS) formulation of 3D gravity [40, 41]. In 
this framework, the EH action is the difference of two CS 
actions; / EH = / C s[A] - / CS [A] with 

I cs [A] = T J TV (A A dA + | A A A A A) (9) 

where A and A are two sZ(2, Revalued connections and 
the CS level k is related to the AdS length £ and Newton’s 
constant Gn through k = . We write asymptotically 

AdS spacetimes in a ‘highest-weight gauge’ [42]. 


A = b 1 (d+a )6 

b = e pL ° 

( 10 a) 

c 

«z = A+ + 

dz = -pL + + ... 

( 10 b) 

Similar expressions hold 

for the barred sector. 

Here 


L±, Lq are the sl( 2 ,R) generators and p is the radial co¬ 
ordinate of AdS. The dots denote terms of lower weight 
which are completely fixed by the EOM. In this gauge, 
the ‘chemical potential’ p corresponds precisely to the 
source of the boundary stress tensor and C is its vev [43] , 
see also [44, 45]. The EOM are the CFT Ward identities: 

— dC = ^9 3 p + 2Ldp + pdC. (11) 

To calculate the 2-point function, we localize the sources 
as in (7) and expand C(z) = C^\z) + eC^{z) + 0(e 2 ), 
where C is the background value, which we take to 
be AdS 3 in the Poincare patch (C^ = L = 0), while 
£(i) corresponds to the linear term in e in ( 8 ). This is 
precisely the 2-point function. From (11) we find 

=Jld 3 ^ 2 \z-z 2 ,z-z 2 ). (12) 

We solve this differential equation using the Green func¬ 
tion of the Laplacean dd in flat space, G{z\ 2l zi 2 ) = 
In (- 212 - 212 )- Then 

£ (1) = ~d 4 zl G(z 12 ) = ^ (13) 

z ~ 12 

gives exactly the 2-point function (4) with central charge 
c = 6 k = 2^7 j compatible with [46]. 

The n-point functions can be obtained similarly. 
We localize the source at n — 1 points, p(z, z) = 
X )”=2 ei5( 2 \z-Zi,z-Zi) := SL 2 On the CFT side, 
this gives the expansion ( T 1 ) M = (T 1 ) + Yli= 2 e i (T 1 T l ) + 
■ ■ • + ( nr =2 d) (T 1 T 2 ... T n ) + 0(e n ). The EOM at order 
e n ~ 1 give 

n 

- dc {n ~ 1) (2£ (n " 2) <% + <5zd£ (n “ 2) ) . (14) 
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Using the Green function G to solve for £ ( " ^ yields 

£ ( " _ 1 ) (*i) = X> ( 4 + — O £ (n_2) (^) • (is) 

i=2 \ Z li Zli / 

The sum of terms in this expression proportional to 
nr = 2 e, is the n-point function (3). All other terms are 
contact terms with two operators localized at the same 
point. We have now derived holographically the 2-point 
function and the CFT recursion relation for the correla¬ 
tors (3) in the CS formulation and therefore know that 
all n-point correlators of the stress tensor match between 
gravity and CFT side. 

The results generalize from the plane to the cylinder 
after taking z = e luJ , with w = it + tp and tp ~ ip + 2n, 
so that (w, ui) ~ (ui + 27r, ui + 2i r). In all correlators 
effectively the only change is that one needs to replace 
the quantities Zij —> 2 sin [(07 — u>j)/2] =: and (3) 

becomes (c,j := cot[(wj — ojj)/ 2 ]) 

(T 1 T 2 ... r)C = E (— + ^ ( T2 • ■ • T ">c (16) 

i =2 Vs li 2 ' 

On the gravity side the relevant vacuum state is now 
global AdS 3 with background values C ^ = |. 


GCFT STRESS TENSOR CORRELATORS 

The asymptotic symmetries of 3D flat space at null 
infinity are governed by the BMS 3 algebra [10, 47] whose 
non-vanishing commutators are (n, m are integers) 

[An; A m ] = (Tl 777 .)A n 4 -m ~f“ l)^n,— m (17a) 

[L n , M m ] = (n- m)M n+m + ^n( n 2 - 1)<5„ _ m (17b) 

This algebra also arises out of non- or ultra-relativistic 
contractions of two copies of the Virasoro algebra (2) 
and is known as 2D Galilean Conformal Algebra (GCA). 
QFTs with GCA symmetries are called GCFTs. If there 
exists a dual formulation of 3D flat space gravity in terms 
of a 2D QFT it has to be a GCFT 2 [10, 12, 13]. 

The stress tensor components Alp, A/p of a GCFT 2 
on the plane can be mode-expanded in terms of L n , M n 
[48, 49], M p = EA/p = En (Ln - (n - 
2)|M„)£ _n_2 . The mapping from the plane to the cylin¬ 
der is given by £ = , 77 = iue lv , with </> ~ + 27 t. On 

the cylinder one has the following mode expansions [49]: 


-M + H = ( 2 M p = E M„e-^ (18a) 

n 

-A + g = C 2 M> - 2£r)Mr = E(^ - inuM n )e~ intp 

n 

(18b) 


For our explorations of flat holography we work on the 
cylinder and consider highest-weight representations of 
the GCFT 2 [50, 51]. We define the vacuum of the theory 
|0) as A„|0) = 0 = M n |0), Vn > —1. Using the above 
relation and its conjugate ((0|L n = 0 = (0| M n , n < 1 ), 
together with (18) and (17), we can calculate any corre¬ 
lation function of the GCFT 2 stress tensor components. 
In particular, the non-vanishing 2-point functions in the 
cylinder representation are given by 

<- M 1 A 2 ) = (A 1 A 2 ) = CL ~ 2 I mTl2 (19) 

Z5i2 ^ s 12 

with the definitions, M l = M(ipi), N l = A(u», <pi), 
Sij '■= 2 sin [(tpi — ipj)/2], Cij := cot [(<p» - <Pj)/2] and 

Tij .— (u^ Uj ) Cij • 

Using the mode expansion and the algebra, the higher- 
point correlators can be derived similarly. Note that (the 
connected part of) any correlator with two or more M- 
insertions vanishes. The first few non-vanishing higher- 
point functions read 


(M 1 A 2 A 3 ) = 2 y 

( 20 a) 

S 12 S 13 S 23 


(A 1 A 2 A 3 )^ Ci r C 2 M E 3 

( 20 b) 

s 12 s 13 s 23 


(M 1 A 2 A 3 A 4 ) — „ 1 CM9 ^ 

( 20 c) 

S 14 S 23 S 12 S 13S24S 3 4 

(A 1 A 2 A 3 A 4 ) - f (7) + CmA4 

( 20 d) 

S 14 S 23 S 12 S 13S24S34 


(M 1 A 2 A 3 A 4 A 5 ) = Ju M ffs( ' 7 ’ ^ 

lll<i<j<5 s fj 

( 20 e) 

(A 1 A 2 A 3 A 4 A 5 ) — 4ci 95 ^ 7 ’ ^ + CM ^ 5 

rii <x< 7<5 s b' 

( 20 f) 

where A 4 = 4 ^( 7)771234 - (ri 234 + n 4 + t 23 )# 4 (7) and 
A 5 = 49 7 g 5 (7, C)^12 3 4 + 49^5(7, C)?72 3 45 — 295(7, C) t 1 2 345- 


In the above, we have defined = Ei <i < j< n T ijj 

e.g. ti 23 = ti 2 +ti 3 + t 23 . The function 34 ( 7 ) [ 55 ( 7 ,C)] is 
again given by (5d) [(5e)], where Zjj has to be replaced by 
in the cross-ratios (5f). The quantity r) pqrs is defined 
as r] pqrs = E'(-l ) 1+i ~ j (ui - Uj)sm(<pk - <pi)/{s 2 pr s 2 qs ), 
and the sum E goes over all six permutations of 
{*, j, k, 1} = 7 r(p, q , r, s ) where i < j and k < l. 


RECURSION RELATIONS 

We derive now GCFT recursion relations involving the 
stress tensor components that encode the corresponding 
Ward identities [51]. The GCFT stress tensor compo¬ 
nents on the cylinder, (18), obey the conservation equa¬ 
tions 

d u (M(<p)(D)= 0 d u (N(u 1 <p)0) = d v (M{<p)0) (21) 
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for any operator (or product of operators) O. The left 
equation (21) looks very similar to the CFT conservation 
equation dz(T(z) O) = 0 [35], which together with the 
similarity of the GCA (17) to the Virasoro algebra (2) 
eventually leads to a recursion relation similar to (16). 



The right equation (21) allows to determine the correla¬ 
tion functions with only IV-insertions in terms of corre¬ 
lations functions with one M-insertion ( 22 ). 


(A r l N 2 ...N n ) = — (M 1 V 2 ...N n ) 

cm 

n 

...N n ) (23) 

i=i 

The second recursion relation (23) can be derived as 
follows. Integrating the right equation (21) and using 
cyclic symmetry of the (M N ... N) correlators we ob¬ 
tain (N 1 ... N n ) = A + £” =1 Uid Vi (M 1 IV 2 ... N n ). The 
quantity A is a w-independent integration constant that 
must be identical to the result in a chiral CFT with cen¬ 
tral charge Cl, since the u- independent part of N contains 
only the Virasoro generators L n . As the (M N ... N) cor¬ 
relators also give chiral CFT results, with central charge 
Cm, we can write A = cl/cm (M 1 N 2 ... N n ). 

We have checked that the recursion relations (22), (23) 
also follow from a suitable contraction of CFT results (16) 
and that they are compatible with all our explicit results 
for correlators, (19)-(20f). Thus we have succeeded in 
deriving the GCFT 2 analog of the CFT recursion relation 
on the cylinder (16). In the remainder of the paper we 
derive the same results holographically. 

FLAT SPACE CORRELATORS 

The 3D flat space EH action can be rewritten as the 
Chern-Simons action (9) for the iso( 2,1) connection A = 
A^M n + A^Ln, where L n , M n (n = 0, ±1) generate the 
global part of the GCA (17) and the trace is defined as 
Tr(A m M n ) = -2^ mn , where 7 = antidiag(l, 1), see 
[31] for our conventions, details and references. Just like 
in the AdS case, we parametrize the radial dependence 
of A by going to a gauge where [27, 31] A = b~ 1 (d+a)b 
with the group element b = exp(^M_) and a = a u du + 
dp dip with it, r, ip the (outgoing) Eddington-Finkelstein 
coordinates. In the absence of chemical potentials the 
connection components take the form a u = M + + Aj M_, 
a v = L + + L- + yp Af- • Here the CS-level k is related 
to Newton’s constant as k = 1/(4 Gjv). The EOM restrict 
the state-dependent functions M, J\f as follows: d u A4 = 
0, d u Af = d v M. These equations are equivalent to the 
GCFT conservation equations (21). 


We proceed as in the AdS case and introduce chem¬ 
ical potentials as sources. The component a v remains 
unaffected, while a u changes to a u + A a u with Aa u = 
— PmAI -|_ — PlL- i- T ciqAIq + cliLq + 02 M_ + a^L_. The 
EOM determine completely the functions cq in terms of 
the chemical potentials //« and pl and restrict the state- 
dependent functions as follows [31]. 

-d u M = kd^pL + pl 9 v M + 2Md v p L (24a) 
-d u M = kd^pM - (1 - PM)d ip M + 2Md v p M 

+ pl dp Af + 2J\fdp> pl (24b) 

Analogous to the AdS case, these equations can be inter¬ 
preted as GCFT Ward identities. 

For our background we choose global Minkowski space, 
A4*' 0 -* = k/2, A A 0 ) =0. To calculate 2-point correlators 
in this Minkowski background we switch on chemical po¬ 
tentials, localized at {u 2: P 2 ), like in the AdS case (7), 
Mm/l = £m/l <5 (2) ( m - u 2 , <p - P2) ■■= £ M /l $■ Here e M/L 
are constants. Plugging these localized chemical poten¬ 
tials into the EOM (24) with A4 = k/2+A4^\ Af = A/’*- 1 ) 
yields the linearized Ward identities. 

d u M (1 ) =- ke L [d^8 + d v S) (25a) 

d u Af (1) = - ke M (d^5 + d v 5) + d v M^ (25b) 

To solve them we use the Green function defined as 
d u d v G(u — u 2 , p — p 2 ) = S. Using the method of images 
yields G(u — u 2 ,tp — ip 2 ) = ln ((u - u 2 ) sin[(</? - ip 2 )/2\) 
and thus we can solve the linearized Ward identities 
(25a), (25b), 

M (i) = 6 k^L = Gk{eM-2e L T V2 ) 

s 12 s 12 

where and are defined below (19). The cl com¬ 
ponent of Ad^ 1 ) (and the cm component of AA 1 )) is 
the 2-point function (M 1 N 2 ) and the cl part of A r( - 1 ) 
is the 2-point function (N 1 N 2 ). They agree precisely 
with the GCFT 2-point functions (19) for cl = 0 and 
cm = 12 /c = 3 /Gjv, compatible with [ 10 , 15]. 

Like in the AdS case we compute the higher-point func¬ 
tions by localizing the sources at multiple points Ui,ipi 
and expanding the vev’s of the deformed theory to order 
e n_1 , which eventually establishes recursive formulas. 

n <~y 

All"- 1 ) = + (27a) 

i—2 VS 1 i 2 ' 

n r 9 

^n-l) = E{(^ + Y dlPi ) ( e M-^ (n_2) + ^ (n ~ 2) ) 
+uM^- + % W"- 2 )) + (27b) 

' S ij S li ^ J 

The cl component of A/"A- 1 ) then corresponds to the 
n-point function ( N 1 N 2 ... N n ). The cm component of 
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j\/-(n-i) C 0 rreS p 0 n( i s to the n-point correlator with one 
M 1 insertion (just as the component of The 

recursion relation (27a) then simply gives (22) with Cm = 
12 k, while the part of (27b) can be simplified to give 

n 

(.N 1 N 2 ... N n ) = Uid Vi (M 1 N 2 ... N n ). (28) 

»=i 

These results show precise agreement with the GCFT 
recursion relations (22), (23), provided that Cl = 0 and 
cm = 12 k = 3/Gjv, again compatible with [ 10 , 15]. 

CONCLUSIONS 

Since the holographic recursion relations (27), (28) co¬ 
incide with the ones that we found in a GCFT (22), (23) 
with the EH values of the central charges, Cl = 0, cm = 
3/Gjv, and the 2 -point correlators coincide as well, we 
conclude that flat space holography is perfectly consis¬ 
tent with the GCFT Ward identities. In other words, all 
holographic n-point correlators of the stress tensor com¬ 
ponents match precisely with the corresponding GCFT 
(connected) correlators. 

It would be interesting to generalize our analy¬ 
sis to other 3D gravity (or gravity-like) theories, see 
e.g. Refs. [19, 20, 40, 52-60] for various intriguing models. 
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